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ABSTRACT 

We present the results from our search for C iv in the intergalactic medium at redshifts z = 5.3 — 6.0. 
We have observed four z ~ 6 QSOs with Keck/NIRSPEC in echelle mode. The data are the most 
sensitive yet taken to search for C iv at these redshifts, being 50% complete at column densities 
log Nq IV ~ 13.4. We find no C iv systems in any of the four sightlines. Taking into account our 
completeness, this translates into a decline in the number density of C IV absorbers in the range 
13.0 < log A^c IV < 15.0 of at least a factor 4.4 (95% confidence) from z ~ 2 — 4.5, where the number 
density is relatively constant. We use our lack of detections to set limits on the slope and normalization 
of the column density distribution at z — 5.3 — 6.0. The rapid evolution of C IV at these redshifts 
suggests that the decrease in the number density may largely be due to ionization effects, in which 
case many of the metals in the z 4.5 IGM could already be in place aX z ^ 5.3, but in a lower 
ionization state. The lack of weak systems in our data, combined with the presence of strong C iv 
absorbers along at least one other sightline, further suggests that there may be large-scale variations 
in the enrichment and/or ionization state of the z ^ 6 IGM, or that C iv absorbers at these redshifts 
are associated with rare, UV-bright star-forming galaxies. 

Subject headings: cosmology: observations — cosmology: early universe — intergalactic medium — 
quasars: absorption lines 



1. INTRODUCTION 

Metal absorption lines are one of our most versatile 
tools for studying the high-redshift Universe. Just as 
the ionization history of the IGM reflects the output of 
ionizing photons from star- forming galaxies and quasars, 
the metal content of the IGM can be used to infer the 
integrated and instantaneous rates of global star forma- 
tion. At very high redshifts, where the direct detection 
of galaxies becomes increasingly challenging, metal ab- 
sorption lines can serve as sign posts to the "typical" 
galaxies that may be responsible for reionizing the Uni- 
verse at z > 6. Perhaps most intriguingly, metal lines can 
also serve as a probe of the reionization process itself. 

A considerable amount of work on the high-redshift 
Universe has focused on cosmic reionization, yet the tim- 
ing and mechanisms by which the IGM becomes per- 
manently ionized remain poorly understood. The three 
main observational constraints on reionization are (1) 
the electron optical depth measured from the polariza- 
tion of the cosmic microwave background (CMB), (2) 
the evolution of the high-redshift Lya forest, and (3) 
the luminosity functions of high-redshift Lya-emitting 
galaxies (LAEs). Five-year IFMj4P measurements of the 
CMB are consistent with an instantaneous reionization 
at Zr 11, but are also compatible with an extended 
reionization process (jPunklev et al.l [2008( 1. The rapid 
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evolution of the mean transmitted Lya flux in the spectra 
of z > 5.7 quasars has been used to inf er that reionization 
may have ended a s recently as z ~ 6 (iBecker et al.ll2001l: 
White et all [20031: iFan et al.ll200l 120061: IGnedin fcFanI 



20061 ). The quasar results remain controversi al, however 
(e.g feongaila 2004; Oh & Furlanctto 2005; Lidz et all 
120061 : [Becker et al.ll2007l ). Lya saturates for even a tiny 
neutral fraction {f^^ ~ 10~^: lFan et al.l[200l ). making it 



extremely challenging to measure the ionization state of 
the IGM at these redshifts directly from the Lya forest. 
Finally, the lack of a sharp decline in the luminosity func- 
tion of LAEs from zsa5.7toz«6.5 suggests that most 
of the volume of the IGM is alre ady h i ghly ionized by z ~ 
6.5. (e.g-lMalho tra & Rhoad^ [200l 120061 : IHu fc Cowid 



l2006l:lDiikstra et al.ii2007t. but see Kashikawa et al. 2006). 
Interpreting the LAE results is somewhat complicated, 
however, as Lya photons may escape through local ly ion- 
ized bubbles (e.g..[Haiman fc Cen 2005; Wvithe fc Loebl 
l2005t iFurlanetto et al.ll2006^ or by redshifting_via galac- 
tic winds before passing through the IGM (ISantosl[2004) . 

Quasar metal absorption lines provide a complimen- 
tary tool for studying the ionization state of the IGM. 
In a reionization scenario where galaxies provide most of 
the ionizing flux, overdense regions of the IGM should be 
the first to become enriched, yet may remain significantly 
neutral until the end of reionization due to the short 
recombination times at high densities. "Forests" of low- 
ionization absorption lines such as O I A1302, Si II A1260, 
and C II A1334 may therefore be present prior to the com- 
pletion of reionization (Oh 2002). Our recent search for 
O I using HIR ES spectra of z > 5 QSOs yielded mixed 
results (Bec ker et al.ll200"6[ ). While an overdensity of low- 
ionization absorbers at z ^ 6 was found along at least one 
sightline, SDSS Jl 148-^5251 (zq = 6.42), no enhance- 
ment was seen towards SDSS J1030+0524 (zq = 6.30), 
which has a complete Gunn-Peterson trough spanning 



2 



Becker et al. 



~ 80 comoving Mpc (| White et al.ll2003f) . These differ- 
ences suggest there may be large-scale variations in the 
enrichment and/or ionization state of the IGM at z ~ 6, 
and perhaps a strong redshift evolution in the ionization 
state of metal absorption systems. The metal-enriched, 
low-ionization regions towards SDSS J1148-f5251 may 
represent the final stages of reionization. Alternatively, 
they may arise by chance from multiple pockets of en- 
riched gas near galaxies aligned along a filament. 

In order to effectively use metal lines as probes 
of reionization, the roles of enrichment and ioniza- 
tion must be separately understood. A quasar sight- 
line that passes through a largely neutral IGM may 
not show an O i forest if the IGM has not yet 
been enriched. At lower redshifts, highly-ionized 
species such as C iv and Si iv have been the pri- 
mary tracers of IG M mctallici ty (e.g., 'Mover fc York! 
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2004HScannapieco et al. 1 120061) . Both the comoving mass 



density of and the column density distribution of C iv 
remain roughly constant over at least 1.5 < z < 4.5 
(|Songailall200l . Since the fraction of carbon in the form 
of C IV varies depending on the ionization state of the 
gas, however, a constant C IV density does not necessarily 
mean that the total IGM metallicity remains unchanged 
(e.g., IQppenheimer fc D ave 2006). Nevertheless, some 
fraction of IGM metals are clearly in place by z ~ 4.5. 

The current constraints on th e abundance of metals 
in the IGM at z > 5 are mixed. iSongailal (|2005D found 
tentative evidence for a decrease in the comoving mass 
density of C IV from z ~ 4 to 5. To measure the metal 
content of the IGM at even redshifts the C iv A1548, 1551 
doublet must be pursued i nto th e ne ar-infrared. Prelim- 
inary searches by iSimcod (|2006f ) and iRyan- Weber et all 
(|2006f ) found evidence for a comoving C IV mass den- 
sity at z ^ 6 that is co nsistent with m easurements at 
z ~ 2 — 4.5. In addition, ISimcod ()200 6^ found tentative 
evidence that the C iv column density distribution may 
also remain constant out to z ^ 6. These works examined 
only two sightlines, however, and were generally sensitive 
only to the rare, high-column density (log A'^c iv ?J 13.8) 
systems that dominate the comoving mass density. A 
complete census of C iv at z 6 requires a broader, 
more sensitive search. 

In this series of papers, we present the results from a 
two-part search for intergalactic metals at z ~ 6. In this 
paper, we describe our survey for C iv performed with 
Keck/NIRSPEC in high-resolution mode. This is the 
most sensitive search for C IV at z > 5.3 yet performed, 
partly due to the high resolution {R « 13,000), which 
is better suited t o search for narrow (6 ~ 10 km s~^, 
iRauch et al.lll996l ) hnes than previous studies. By cover- 
ing multiple sightlines with high sensitivity, we are able 
to place strong constraints on the abundance of the kinds 
of highly-ionized absorbers that are common at lower 
redshifts. In the following paper, we will expand our 
search for Iqw-ioni zation systems initially presented in 
iBecker et al.l ()2OO60 . The combined results will provide 
unique constrains not only on the ionization state of the 
IGM, but on the star formation history of the Universe 
at z > 5. 



The rest of the paper is organized as follows. The data 
are presented in §2. In §3 we give the results of our 
search for C IV, along with our completeness estimates. 
We show that the number density of C iv absorbers at 
z > 5.3 has declined sharply from z ~ 2 — 4.5 in §4, 
and give constraints on the column density distribution 
and the contribution of C iv to the closure density. In 
§5 we discuss possible reasons for the decrease in C iv 
with redshift and the significance of our lack of detections 
in light of the strong systems previously reported along 
other sightlines. Finally, our conclusions are presented 
in §6. Throughout this paper, we assu me Sim = 0.274, 
p,\ = 0.726, and Hq — 70.5 km s~^ (|Komatsu et al.l 
[200l) . 

2. THE DATA 
2.1. NIRSPEC Echelle spectra 

We obtained high-resolution, ne ar-IR spectra of four 
z ~ 6 QSOs with Keck/NIRSPEC (|McLean et al.|[l998l) 
between February 2007 and February 2008. The spec- 
trograph was used in echelle mode in the NIRSPEC- 1 
(Y) band, which provides continuous wavelength cover- 
age over A « 0.95 — 1.13 /im in a single setting. Our 
observations are summarized in Table [TJ Individual ex- 
posures were 1800 seconds, with total integration times 
of 11-12 hours per object. 

Data reduction was performed using a custom set of 
IDL routines tailored to the technical challenges of ob- 
serving faint targets in the infrared. The NIRSPEC de- 
tector, a 1024x1024 ALLADIN-3 InSb array, has good 
sensitivity in the near infrared but also a high dark cur- 
rent. The sky emission lines are widely separated in 
echelle mode, leaving the dark current as the dominant 
source of noise provided that other sources such as vari- 
able bias levels, hot pixels, and numerous cosmetic blem- 
ishes can be properly removed. 

Rather than nod-subtract our exposures, which auto- 
matically multiplies the noise from the sky and the dark 
current by a factor of root two, we compiled a low noise 
dark frame by taking a large set of dark exposures with 
the same exposure time as our science frames. This was 
used to remove dark current features and other blemishes 
prior to sky subtraction. The exposures were flat fielded 
using internal flats. In order to correct for pixel-to-pixel 
variations separately from the high-frequency fringing in- 
troduced by interference between the order sorting and 
blocking filters, the fringe pattern was modeled and re- 
moved from the two dimensional flats. Each science ex- 
posure was then divided by the "pixel" flat, as well as 
by the fringe model after it was scaled in intensity and 
shifted on the detector to match the exposure. 

The sky was modeled using a b-spline fit and removed 
from ea ch exposure using optimal sky subtraction tech- 
niques (jKelsonl 120031) . The residual bias level from each 
of the thirty-two independent readout channels was mod- 
eled and removed both at the start of the reduction and 
again after a first-pass sky subtraction. Similar care was 
taken to remove the residual bias in the dark frames. 

In order to maximize the signal-to-noise ratios of the 
final spectra, we performed simultaneous multiple spec- 
tral extraction on all exposures for a given object. The 
object spatial profile was modeled independently for each 
exposure. A response function derived from a standard 
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TABLE 1 

Summary of NIRSPEC Echelle Observations 



QSO=' rcdshift tcxp (hrs) S/N"" Azc iv" AXq iv'^ 

SDSS J0002+2550 5.82 12.0 11.4- 16.7 0.458 2.33 

SDSS J0818+1722 6.00 11.0 10.3- 14.2 0.638 3.29 

SDSS J0836+0054 5.80 11.5 14.5 - 21.5 0.438 2.23 

SDSS J1148+5251 6.42 12.0 10.9- 15.6 0.679 3.44 

^ QSO names given as SDSS Jhhmm+ddmm 

^ Middle 50% range of signal-to-noise per resolution element values for pixels 

within the wavelength interval used to search for C IV 
Total redshift interval used to search for C IV 
Total absorption path length interval used to search for C IV 



star exposure was also applied to the counts in the 2-D 
frames. A single 1-D spectrum was then extracted from 
all exposures of a given obje ct simultaneously using op- 
timal weighting ()Hornelll986f ). This technique allowed us 
to best reject spurious pixels while preserving as much 
information as possible from the 2-D arrays, which were 
individually very noisy even after careful reduction. The 
1-D spectra were then normalized using a slowly-varying 
spline fit to the continuum. Telluric absorption correc- 
tion was done using a model of the atmosp heric trans- 
mission in the Y band (|Hinkle et al.l 120031 ). This pro- 
duced better results than using our standard star spec- 
tra, which contained some residual high-frequency ring- 
ing as a result of the fringing. Even so, we have restricted 
our analysis to the center of the transmission window 
(9850 - 10920 A), where there is very little atmospheric 
absorption. Finally, we scaled the formal error arrays so 
that the distribution of deviations of the flux from the 
continuum fit, normalized by the error array, was well fit 
by a Gaussian with unit variance for each object. The 
typical resulting increase in the error array was ~ 30%, 
which is a reasonable estimate for the amount of noise 
that could not be well characterized in the 2-D frames. 

Representative sections of the final spectra are shown 
in Figure [TJ The data tend to become noisy at the ends 
of the echelle orders, which are marked by a step in 
the error array. The measured resolution is i? w 13000 
(FWHM « 23 km s~^). Narrow metal lines can clearly 
be seen, though these are at much lower redshifts. Signif- 
icantly, the noise properties of the data are well-behaved, 
with a scatter as a function of wavelength that is well- 
described by the scaled error array. This is essential for 
correctly estimating our sensitivity. 

2.2. Optical spectra 

We have also acquired high-resolution optical spectra 
of z ^ 6 Q SOs. Some of these data were presented in 
iBecker et a l. (2006, 2007), and the new observations will 
be discussed more thoroughly in the r iext p aper. Here 
we win use Keck HIRES ( Vogt et all 11994 ) spectra of 
SDSS J0836-h0524 and SDSS J1030-h0524 [zq = 6.30) 
to help evaluate candidate C iv detections along those 
sightlines. The data were reduced using a custom IDL 
pipeline featuring optimal sky subtraction. As with the 
NIRSPEC data, simultaneous optimal extraction was 
also used to maximize the signal-to-noise ratio of the fi- 
nal spectra. Prior to extraction, the 2-D counts were 
flux calibrated using a response function generated from 
a standard star. Telluric absorption correction using a 
standard star was also applied to the 2D frames. Veloc- 



ity offsets between the standard and the science frames, 
as well as variations in the strength of the H2O and O2 
absorption bands were determined from preliminary ex- 
tractions of groups of individual exposures taken near 
the same time. Since SDSS J1030+0524 was observed on 
multiple runs at different times of the year, shifts in the 
relative velocity of the object allowed residuals from at- 
mospheric absorption and emission to be especially well 
rejected. 

3. RESULTS OF C iv SEARCH 

We searched for C iv doublets in our NIRSPEC data 
both by eye and using an automated technique described 
below. Although the spectra were corrected for telluric 
absorption, we limited our search to the regions between 
9820 A and 10920 A (5.34 < zc„ < 6.04), which is 
largely free of atmospheric absorption. The total redshift 
interval along each sightline is given in Table [TJ In cases 
where the QSO C iv emission line fell in the bandpass, 
we included the entire region blueward of the line center. 
While it is customary to exclude a narrow region near the 
QSO redshift in order to avoid proximity effects, this re- 
gion is a lso typical ly enhanced with C iv absorbers (e.g., 
lVestergaardli2003i: iNestor et al.ll2008[ ). Our conclusions 
below are therefore most likely enhanced by including it 
in our analysis. In any case, the proximity region com- 
prises only a small part (~ 10%) of the total path length 
towards two of our four QSOs. 

Both the manual search and the follow-up automated 
search revealed only a single candidate C iv absorber, 
which is at z = 5.4091 towards SDSS J0836-f0054. 
This system is plotted in Figure [2] Voigt profile fit- 
ting produced a best-fit column density log = 13.06 ± 
0.12 (cm~^) and a Doppler parameter b — 5.8 ± 
6.3 km s~^(i.e., unresolved). Using the HIRES spectrum 
of SDSS J0836-I-0054, we can measure Si IV and H i Lya 
and Ly/3 at the same redshift. There is no detectable 
Si IV, although this does not necessarily rule out a le- 
gitimate C IV absorber. More significantly, the HIRES 
data show transmission in Lya and Ly/3 at this redshift. 
Since C iv absorbers of this column density are almost al- 
ways accompanied by s aturated Lyg absorption, at least 
at low er redshifts (e.g.. lElhson et al]ll999l : TSimcoe et al.l 
I2OO4 but see Schaye et al. 2007), we can safely conclude 
that this system is not real. Thus, both our manual and 
automated searches produced no C IV detections. 

3.1. Completeness 

To assess the significance of our lack of C iv detections 
we have estimated our completeness as a function of col- 
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Fig. 1. — Example sections of the final 1-D normalized NIRSPEC echelle spectra. The formal error arrays are shown in red. Genuine 
absorption lines are identified with tick marks. The lines in the spectrum of SDSS J0002+2550 are Fe II A2344, 2374 at z = 3.404. Those in 
the spectrum of SDSS J1148+5251 are three sets of Mg II A2796,2804 doublets at z 2.77. Features in the fiux that coincide with spikes 
in the error arrays are typically residuals from the subtraction of strong skylines. 



umn density using a Monte Carlo method. Artificial C iv 
doublets were inserted randomly into the data and an au- 
tomated algorithm was used to determine whether these 
systems could be detected. We generated 10* doublets 
for each sightline, ranging in column density from 10^^ 
to 10^^ cm^^. Doppler par ameters were select ed from 
the measured distribution of I Ranch et al.l (|1996l ). 

In order to be considered a detection, we first re- 
quired the flux decrement to exceed the local error at 
the wavelength of both lines of the doublet after the 
spectra were smoothed by the instrumental resolution 
(FWHM = 23 km s~i). We then fit independent Gaus- 
sians to each line (using the unsmoothed data), and eval- 
uated the reliability of the fits and the agreement be- 
tween the two lines. Our tolerances were stricter for the 
stronger A1548 transition. For a satisfactory fit, the ab- 
solute velocity uncertainty was required to be less than 
0.5 (1.0) times the resolution for the A1548 (A1551) tran- 
sition. We further required that the fitted flux decrement 
be greater than 3.0 (2.0) times the fit uncertainty, and the 
line width be greater than both 3.0 (1.0) times the fit un- 
certainty and 0.5 times the resolution. A "detection" was 
only counted when the line fits were also consistent with 



being a doublet. In order to exclude potential false posi- 
tives, we set stricter requirements for weaker lines, since 
strong lines are typically easy to identify by eye. When 
the peak optical depth of the A1548 transition was less 
(greater) than 1.0, we required that the velocity offset 
between the two fits be less than both 1.0 (2.0) times the 
resolution and 2.0 (3.0) time the combined uncertainty 
in the fitted velocities. The line widths were required to 
agree to within half the width of the narrower line, or 
three times the combined uncertainty. In all cases, how- 
ever, the line widths could not disagree by more than a 
factor of two. Finally, if the peak optical depth of one or 
both lines was less than one, we required the ratio of the 

optical depths ■'■ai548/''a'i551 between 0.8 and 3.0. 

Defining criteria for an automated search is a some- 
what subjective process, as selection thresholds and 
goodness-of-fit requirements will vary depending on the 
line fitting method. The above procedure was designed to 
give the same result as a visual evaluation of the synthetic 
lines as often as possible. In general, this meant apply- 
ing greater scrutiny to weaker systems in order to avoid 
false positives. We found very good agreement overall 
between the fraction of systems identified by eye and by 
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Fig. 2. — Candidate C IV system at 2 = 5.409 towards 
SDSS J0836+0054. Histograms show the normalized flux, which 
is from NIRSPEC for C IV and HIRES (binned to 4.2 km s^^ for 
display) for Si IV and H I. Red lines are the formal error arrays. 
The blue line is a Voigt profile fit to the candidate C IV absorp- 
tion, with A^c IV = 10 cm~^. The presence of H I Lyo and 
Ly/3 transmission at the same redshift strongly indicates that this 
is a false detection. 



the automated method. 

Our completeness results are shown in Figure [31 We 
are most complete at logiVc iv > 13.7, though we are still 
50% complete in the overall sample at log A^c iv ~ 13.4, 
and even retain some sensitivity at logiVc iv = 13.0. As 
discussed below, the fact that we are significantly com- 
plete over a wide range of column densities allows us to 
put strong constraints on the column density distribution 
of C IV absorbers. 

We also evaluated the number of false positive de- 
tections we should expect. Using the same automated 
search method, we searched the data for "fake" doublets 
with velocity separations that were offset from the C iv 
doublet separation. In each trial we adjusted the doublet 
separation by twice the instrumental resolution, omitting 
separations that were similar to true doublets such as 
Mg II A2796, 2804, or other plausible combinations of in- 
tervening absorption lines. In twenty trials we detected 
seven false doublets, which gives an expectation value of 
0.35 for the entire data set. This is consistent with the 
single false C iv detection described above, although in 
that case we were able to rule it out based on the presence 
of transmitted flux in the Lya forest. 



4. CONSTRAINTS ON C iv IN THE Z ~ 6 IGM 

4.1. Change in the number density of absorbers from 

z ~ 3 

We can use our completeness estimates to calculate 
the number of detections we should expect if the column 
density distribution has not evolved from lower redshifts. 
It is conventional to fit a power law of the form 
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Here, X is the absorption path length interval, de- 
fined such that, for a non-evolving population of sources, 
dAf/dX remains constant with redshift. For = 1, 



dX 
dz 



[l + zf 



V^Mi 
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Fitting a large sample of absorbers at z = 2.90 — 3.54, 
iSongailal (|200lD found a slope a = 1.8 ± 0.1 for sys- 
tems with logiVciv > 13.0, with a normalization B = 
j^Q- 12.67 (corrected for our adopted cosmology) at Nq — 
lO^'^ *' cm~2. Indeed, they found that this distribution 
is nea rly invariant out to at least z k, 4.5. ISongailal 
(|2005l ) found a slightly shallower slope (a = 1.7) and 
a higher normalization when fitting the distribution of 
"pseudoclouds" identified by an automated method. We 
will compare our results to the[Songaila (2001) fit, how- 
ever, since the search method is more similar to the one 
used here. 

Figure |4] shows our mean completeness func tion con- 
volved with the column density distribution of ISongailal 
()2001f ). We would expect the largest number of detec- 
tions near log A^c iv ~ 13.5, but there are significant tails 
towards both ends of the column density range. In total, 
we would expect to detect 13.2 systems with log iVc iv be- 
tween 13.0 and 15.0 along all four sightlines if f{Nc iv) 
was unchanged from z ~ 3. For zero detections, the 95% 
and 99% 1-sided confidence limits on the mean expected 
number are 3.0 and 4.6, respectively (e.g., Gehrels 1986). 
This translates into a 95% confidence limit of at least a 
factor of 4.4 decline in the number density of C iv sys- 
tems from z 3 to z = 5.3 — 6.0, or a 99% Hmit o f 
a factor of 2.9. We note that with the ISongailal ()2005l ) 
f{Nc iv) we would have expected 21 systems. Similarly, 
using f{Nc IV ) from a fit to the C iv distribut ion in 
nine high-quality HIRES spectra from Boksenbcr g et al.l 
(2003), w e wou ld have expected 17 lines. Adopting the 
iSongailal (j2001h column density distribution may there- 
fore be somewhat conservative. 

4.2. Previously studied sightlines 

Given the large drop in the number density of ab- 
sorbers we observe, it is worth reviewing t he resul t s frorn 
the initial searches for C iv at 2 > 5.3. iSimcod ()2006f ) 
obtained Gemini/GNIRS spectra of two z ~ 6 QSOs, 
SDSS J1306+0 356 (zo = 6.00) a n d SP SS J1030+0524 
(zQ = 6.27). iRyan- Weber et all (|2006l ) also observed 
SDSS J 1030-1-0524 with VLT/ISAAC. Simcoe found a 
number of low-column density absorbers along both 
sightlines. As the author points out, however, the ma- 
jority of these detections are near the detection limit and 
are consistent with the expected number of false detec- 
tions. They found at least one secure absorber, which 
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Fig. 3. — Percentage of recovered artificial C IV doublets as a function of column density. Histograms give the completeness for individual 
sightlines. The dotted line is the path length-weighted mean for the entire sample. The completeness does not go to 100%, even for large 
column densities, because of the non-negligible chance that of one of the C IV components may coincide with a lower-redshift absorption 
line, a residual from a subtracted skyline, or some other unusable part of the spectrum. 
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Fig. 4. — Expected C IV detection rate for a column density dis- 
tribution that is unchanged from z 2 — 4.5. We computed the 
expected number of C IV systems per unit redshift pat h interval 
per lo g column density bin by integrating f(Nc iv) from [5oneailal 
11200 II) against our path length-weighted mean completeness (Fig- 
ure [SJ. Over our entire sample, we would expect to detect 13.2 
systems with 13.0 < logA^c iv < 15.0. 



in the GNIRS spectrum). They further detect a system 
at z = 5.724 that is unclear in the GNIRS data. We 
have checked the authenticity of the z = 5.724 system 
by searching for related lines in our HIRES spectrum of 
SDSS J1030-H0524. The optical data, shown in FigureEl 
show Si IV absorption at the redshift of the candidate 
C IV. We conf irm the large 60 km s~ ^) velocity width 
measured by IRvan- Weber et al.l ()2006D . At HIRES res- 
olution, the Si IV absorption is clearly asymmetric and 
best fit using multiple components (see Table 2). As- 
suming that the C IV absorption has the same kinematic 
structure and relative component strength as the Si iv, 
we were able to reproduce the Ryan-Weber et al. fit 
using a total C iv column density very similar to their 
measured log A'civ — 14: A. Thus, there do not seem to be 

any hidden saturation effects in the ISAAC data^^ 

In s umm ary, in the two sig h tlines studied by iSimcod 
(|2006f) and IRvan- Weber et al.l (|2006f ) there are two con- 
firmed z ~ 6 C IV systems. These lie along a single 
sightline, and they are both strong (log A^n tv > 13 - 8) sys - 
tems. The weaker candidates reported bv ISimcod ()2006l ) 
are broadly consistent with their expected false positive 
rate. 



occurs towards SDSS J1030+0524 at z = 5.829 with 
logN^l^^ = 13.8. Ryan- Weber et al. also detected this 
system, albeit marginally (they measure a much higher 
column density, but this may be due to blending of the 
A1551 transition with an unrelated line that is visible 



4.3. Constraints on the column density distribution 

We can set joint constrains on the slope and normaliza- 
tion of a power-law distribution of C IV column densities. 
For a given slope a, we have calculated the 95% limits 
on B by adjusting the normalization such that we would 
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Fig. 5. — HIRES data showing the Si IV and H I absorption 
profiles for the can didate C IV system at z = 5.724 towards 
SDSS J1030+0524 IIHvan- Weber et al.|[200g) . The spectra have 
been binned using 4.2 km s~^ pixels for display. The blue line is 
a Voigt profile fit to the Si IV using three components (Table [2]l . 
While the data are too noisy to determine the exact kinematic 
structure, the Si IV is clearly not saturated, and therefore the to- 
tal fitted column density does not change significantly when more 
components are added. 



TABLE 2 

VoiGT Profile Fit to Si iv at z = 5.724 towards 
SDSS ,11030+0524 



z 


b (km 


3 




log A'' (cm' 




5.7235623 ± 0.000049 


11.8 


± 


3.6 


12.97 


± 


0.10 


5.7242061 ± 0.000063 


9.2 


± 


4.2 


12.98 


± 


0.17 


5.7249523 ± 0.000042 


16.2 


± 


2.7 


13.47 


± 


0.06 



Note. — These are the best fitting parameters for 
a three component fit to Si IV in the HIRES data. 
Increasing the number of components does not sig- 
nificantly change the total column density. 



expect to detect 3.0 C iv systems along our four sight- 
lines. In calculating i?, we take into account our mea- 
sured completeness, and restrict ourselves to the number 
of expected systems between logiVc iv = 13.0 and 15.0, 
for reasons discussed below. The results are shown in 
Figure [6l As stated above, if a remains unchanged from 
z 3, then the normalization must decrease by at least 
a factor of ~ 4. Alternatively, iv) could steepen 

(larger a) while maintaining a somewhat higher normal- 
ization. This would effectively "hide" a greater fraction 
of the expected systems below our detection limit. It is 
also possible that j{Nc iv) could become more shallow, 
shifting the distribution towards stronger lines, provided 
that the normalization decreases more sharply. In Fig- 
ure[6]we have only shown the results for a > 1. Formally, 
a could be less than one, but this would imply more ab- 
sorbers per decade in iv with increasing Nq which 
is not observed for any other class of absorbers. We fur- 
ther note that low values of a strongly require the as- 
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Fig. 6. — Limits on the column density distribution, f{Nc iv) = 
BNg"^, of C IV absorbers a,t z = 5.3 - 6.0 in the range 13.0 < 
log A^c IV < 15.0. The solid shaded regions show the range of nor- 
malizations allowed at the 95% confidence level for a given power- 
law slope based on the lack of detections in our NIRSPEC data. 
Separately, the line-filled regions are excluded at the 95% level 
based on the fact that there are at least two known systems with 
log Ac IV ^ 13.8 along six lines of sight (see text). Symbols show 
fitted parameters for /(Ac i v) at z ~ 3, corr ected for our ad opted 
cosmology; square: Songaila 11200 II) . triangle: ISoneailal II2005I) . cir- 
cle: B oksenberg et al.i (i2003i ). 

sumed high-end cutoff in A^c iv at ip is.o cm~ ^. 

We can further u s e the detections of lSimcoj (|2006D and 
iRvan-Weber et al] ()2006l ) to set lower limits on the nor- 
malization of f{Nc IV )■ The total path length surveyed 
by those authors and in this work is AX = 17.4. For two 
detections, the 1-sided 95% lower limit on the expected 
mean number is 0.36. In order to set a lower limit on 
B, we assume 100% completeness and adjust the nor- 
malization such that we would expect at least this many 
systems with log A'c iv = 13.8 — 15.0 along all six sight- 
lines. The results are shown in Figure [B] Including this 
constraint still leaves a significant degeneracy between a 
and B, but it greatly reduces the range of allowed nor- 
malizations for a given power-law slope. 

4.4. Limits on flc ,y 

Our limits on the total C iv mass density are nec- 
essarily weaker than our constraints on /(A^civ)- For 
a < 2, the integrated mass density will be dominated 
by the highest column density absorbers, which may be 
too rare to appear in our data. Nevertheless, we can set 
limits on Qq iv by integrating over the allowed column 
density distributions. The contribution of C IV to the 
closure density pc can be calculated as 

nc IV = — / A^c IV f{Nc iv) dNc IV , (3) 

Pc C Jjv„i„ 

where mc iv is the mass of the C iv ion. The choice 
of integration limits is critical, but there are natural 
choices for C IV. At lower redshifts, f{Nc iv) tends to 
become signifi cantly shallower at logiVciv ^ 13.0 (e.g., 
ISoneailal 120051 ) . Ii itegrating the c ompleteness-corrected 
fits to /(A^civ) of lSongailar()2005f ) yields a contribution 
to r^c IV below 10^^ ° cm~2 that is < 10% of the con- 
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Fig. 7. — Limits the contribution of C IV to the closure density 
a.t z = 5.3 — 6.0, calculated by integrating over the allowed column 
density distributions in Figure|6] The solid shaded region shows the 
range of Qq iv values allowed at the 95% confidence level for a given 
slope a, based on the lack of detections in our NIRSPEC data. Sep- 
arately, the line-filled regions are excluded at the 95% level based 
on the fact that there are at least two known log A^c iv ^ 13.8 sys- 
tems along six lines of sight. The values for Qq jv were calculated 
by integrating Eq.[3]over lO^"^" cm~^ < A^c iv < lO^^" cm~^ (see 
text). The horizontal dotted line shows Q q iv at z '--^ 2 — 4.5 cal- 
culated by integrating over the iSongaila (20011) f{Nc iv) between 
the same bounds, which produces good agreement with published 
observations. 

tribution between 10^^ ° and 10^^ '^ cm^^. At the other 
end, systems with logiVc iv > 15.0 are not typica l ly de- 
tected, even in large surveys fe.g.. lSongailall200l I2005L 
but see Scannapieco et al. 2006). There is likely to be 
another break in the column density distribution above 
log 7V"n IV ~ 15.0, at least at moderate redshifts. Songaila 
(|2001l ) would have expected to detect ~ 5 systems with 
15 < logiVciv < 16 between z = 2.0 and z = 4.5 if 
their power-law distribution held to these column densi- 
ties. At z ~ 2 — 4.5, the measured fie t v ^ 3 x 1 0~^ is 
well-reproduced by integrating over the ISongailal (IMlh 
f{Nciv) between 10" ° and lO^^" cm'^. We therefore 

adopt TVmin = 10^3.0 ^^-2 g^j^j j^^^^ ^ ;^Ql5.0 

our integration limits. 

The results of integrating over /(iVc iv) using the al- 
lowed values of the power-law slope and normalization 
are shown in Figure [71 The corresponding value at 
z ^ 2 — 4 5, ob tained by integrating over f{Nc iv) from 
ISongailal (|200H ). is shown as a horizontal line. Our data 
at least demonstrate that fie iv(^ = 5.3 — 6.0) must be 
< 1.3 X 10-^ or < 0.5 nc iv(z ~ 2 - 4.5). As expected, 
the largest allowed values of fie iv are for very shallow 
/{Nq iv), which permit there to be rare, strong systems 
without also requiring a large number of weaker systems. 
If a remains unchanged from z ^ 2 — 4.5, however, then 
flc IV is limited to < 6.0 x 10^^, or a factor of ~ 4 less 
than ilc iv{z ^2 — 4.5). Separately, Sic iv(-z = 5.3 — 6.0) 
must be > 1 X 10^^ in light of the two detected systems 
towards SDSS J1030-|~0524. These limits are compared 
to measurements from the literature in Figure[51 We note 
that our limits on fie iv based on the allowed shape of 
the column density distribution are consistent with the 
results from a broader search for strong C IV systems 



et al. 



10.0 



X 




2 3 4 5 6 

z 

Fig. 8. — IV as a function of redshift. The downward-pointing 
arrow shows the 1-sided 95% upper limit on f!c tm sk, z = 5.3 — 6.0 
set by the lack of detections in our NIRSPEC data (see Fig- 
ure (Tjl. The limit assumes a power-law column density distribu- 
tion, /(Afc iv) oc ^c%, with a > 1.0 for 13.0 < logiVc iv < 15.0. 
Separately, the upward-pointing arrow shows the 95% lower limit 
on Q.Q IV for a > 1.0 set by the fact that there are at least two 
known log A^c iv > 13.8 systems along six lines of sight. Symbols 
are measurements from the literature, corrected for our adopted 
cosmology, with 2-sided 90% vertical error bars; small squares: 
ISongaila (2001), filled triangles: iSon gaila (20 03), o pen c ircles: 
[Boksenbere et al.l (120031) . open triangle: iPettini et al.l 120031 '). 

along a larger number of z ~ 6 sightlines (Ryan- Weber 
et al., in prep). 

5. DISCUSSION 

Our primary goal in this paper has been to measure 
the evolution of C IV absorbers aX z — 5.3 — 6.0. We 
have shown that the weak (logiVc iv ^ 13.0) absorption 
systems that are ubiquitous at z < 4.5 are less abundant 
by at least a factor of four at z > 5.3. At face value, 
this might suggest that the metal content of the IGM is 
changing rapidly at these redshifts. Several factors, how- 
ever, must be considered when translating the observed 
incidence of C IV systems into a global IGM metallicity. 
The number density of highly ionized metal absorption 
systems at a given redshift will depend not only on the 
mean metallicity of the IGM, but also on the the spatial 
distribution of the metals and on the ionization state of 
the enriched gas. We will have more leverage to assess 
these factors in a subsequent paper, in which we exam- 
ine the role of low-ionization metal absorption systems at 
z ~ 6 (Becker et al., in prep). For now, however, we can 
offer some speculation based on the C iv results alone. 

Perhaps the most striking feature of the high-redshift 
C IV evolution is how rapidly it proceeds at z > 4.5. As 
noted above, both the integrated comoving mass density 
(fic iv) and the column density distribution (/(A^civ)) 
remain nearly invariant over 2 < z < 4.5. Our redshift 
coverage extends down to z = 5.3, which means that 
much of the buildup of C iv, at least for the weaker 
systems, must happen over an interval Az < 1, or 
Ai < 300 Myr. By comparison, neither the comoving 
star formation rate density nor the stellar mass density 
is likely to increase by more than a factor of three from 
z ~ 5.3 - 6 to z ^ 4.5 (Bouwons et al. 2()ol lYan et al.l 
[2006L IStark et al.l[2007t l Evlcs ct al. 2001). 

lOppenheimer fc Dave! (2006, 2008) have recently ad- 
dressed the observed C iv evolution using numerical 
simulations of large scale structure, including various 
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prescriptions for feedback from galactic winds. They 
find that a relatively constant flc iv and f{Nc iv) over 
2: ~ 2 — 5 can be produced by a steadily increasing IGM 
metallicity offset with a declining C iv/C ratio due to en- 
ergy input from winds. Significantly, they find that in or- 
der to produce enough systems a,t z = 3.0 — 4.5, the met- 
als must largely be in place by z = 4.5 — 6.0. Their best- 
fitting model ('vzw' from iQppenheimer fc Pavel 120061 ) 
uses momentum-driven winds whose strength is tuned 
to strike the best balance between metal enrichment and 
energy injection. Convolving their f{Nc iv) from this 
model at z 4.5—6.0 with our completeness, and making 
a first order correction based on the fact that their Qq iv 
declines by roughly a factor of two from z ^ 4.5 — 5.3 to 
z ~ 5.3—6.0, we would expect to see ~ 13 C iv systems in 
our data. This is strongly excluded, which suggests that 
the mean metallicity, the spatial distribution of metals, 
and/or the ionization fraction of C iv at z > 5.3 must be 
evolving more rapidly than in the ir simulations. 

An important feature of the IQppenheimer &: Pavel 
(|2006f ) results is that C IV traces incre asingly smalle r 
overdensities at higher redshifts (see also iSimcod 120061 ) . 
The observed dropoff in C IV may therefore indicate that 
most of the metals at z > 5.3 are confined to over- 
densities greater than those where C IV is a sensitive 
tracer of metallicity (—1 5 , fogp/p ^ 1 at 2; ~ 5.5; 
IQppenheimer fc Pavel 120061 ). The metallicity gradient 
might have to be fairly steep in order to suddenly recover 
the observed number density of C IV systems at z ~ 4.5 
by shifting the C iv "window" to larger overdensities. If 
the metal-enriched regions are marginally self-shielded at 
z > 5.3, however, then a moderate increase in the UV 
background may be sufficient to produce a large increase 
in the number of C iv systems over a relatively small 
Az. A hardening of the UV backgrou nd could also pro- 
duce a n increase in the C iv/C ratio. iMadau &: Haardtl 
(|2008f ) have recently pointed out that absorption in the 
Lyman series of He 11 in the IGM can substantially de- 
crease the intensity of the UV background at energies 
near the ionization potential of C ill (3.5 Ryd). If the 
volume filling fraction of He 11 is significantly smaller at 
z ^ 4.5 than at z > 5.3, as would be ex pected for an 
extended He 11 reionization process (e.g., iBolton et al.l 
120081: iMcQuinn et~aLll2008D . then the shape of the UV 
background may be more favorable to C iv at z ~ 4.5 
than at z > 5.3. 

Given t he strong C IV systems i dentified by ISimcod 
((2006) and iRvan- Weber et al.l (|2006l ). it is especially sur- 
prising that we did not detect any weak systems. The two 
C IV systems towards SPSS J1030-I-0524 are marginally 
consistent with a power-law column density distribution 
under the constraints provided by our data. For example, 
for our 95% upper limit on the normalization of /(A^c iv) 
for a = 1.8 (Figured]), there is a 6% chance that a sin- 
gle sightline spanning AX = 3.1 would contain at least 
two systems with 13.8 < logA"civ < 15.0. Among the 
six sightlines included here and in the previous studies, 
however, there would be a 30% chance that at least one 
sightline would contain two or more strong systems. This 
increases to 63% for a = 1.0, but in either case the lack 
of detections along our four sightines would have been 
highly unlucky. It is intriguing, therefore, to consider 
the possibility that there may be large-scale enrichment 
and/or ionization fiuctuations in the z ~ 6 IGM. 



As noted by ISimcod ()2006D . there is good evi- 
dence that strong C IV systems at lower redshifts are 
closel y associa ted with actively star-forming galaxies 
( Ade lberger et al. 2003; Simcoe et al. 2006). If weak C iv 
systems are absent from the z 6 IGM because most of 
the metal-enriched gas is in a lower ionization state, then 
the rare strong C IV systems may point to regions near 
vigorously star-forming galaxies, where the metallicity 
and the local ionizing background are enhanced. There 
is tentative evidence that this is the case from a search 
for z ~ 6 Lyman break galaxies carried out by Kim et al] 
(2OO8). They used HST/ACS to search for i-dropouts in 
the fields of five z ~ 6 QSQ. Compared to the GQQPS 
field, they find overdensities in only two fields, with the 
strongest overdensity towards SPSS J1030-I-0524 (first 
reported by Stiavelli et al. 2005). The i-dropouts in that 
field appear spatially clustered, and their i — z colors 
are consistent with star-forming galaxies at z ~ 5.8 
(|Bouwens et al.ll2006() . This suggests that the two strong 
C IV absorbers seen towards SPSS J1030-I-0524 may be 
associated with outflows from UV-bright galaxies in a 
foreground filament. 

In contrast, an underdensity of i-dropouts is seen to- 
wards SPSS J1148+5251, whose sightlines does not con- 
tain any z > 5.3 C iv absorbers. This sight hne is known 
to co ntain multiple z '--^ 6 Q i systems (jBecker et al.l 
I200I) . The presence of metals in a low-ionization state, 
combined with the lack of UV-bright galaxies, suggests 
that this sightline was enriched at z > 6, but lacks 
a sufficient amount of ongoing star formation to keep 
the enriched gas highly ionized. We will return to this 
point in a subsequent paper. Finally, we note that the 
SPSS J08 36-I-0054 field co ntains an overdensity of i- 
dropouts ( Zheng et al. 2006) yet does not show any C iv 
or Q I (Becker ct al. 2006). If these galaxies are embed- 
ded in a protocluster surrounding the radio-loud QSQ, 
then hard UV photons from the QSQ itself may be over- 
ionizing any metal-enriched gas such that it does not ap- 
pear as C IV. Spectroscopic follow-up of the z-dropouts 
in QSQ fields will help determine whether these scenarios 
are correct. 

6. CONCLUSIONS 

We have presented a survey for C iv at z = 5.3 — 
6.0 towards four z ~ 6 QSQs using Keck/NIRSPEC in 
echelle mode. This is the most sensitive search for C iv at 
these redshifts to date, yet we find no C iv systems in any 
of our spectra. Based on our com pleteness estima tes and 
the column density distribution of lSongailal ()2001l ) , which 
remains largely invariant over 2 < z < 4.5, we would 
have expected to detect roughly 13 systems in the range 
13.0 < logA'^civ < 14.5 in the case of no evolution out 
to z '--^ 6. Qur lack of detections implies that the number 
density of C iv absorbers to which our data are sensitive 
declines by at least a factor of w 4.4 (95% confidence) 
from z ~ 4.5. 

Qur null result places strong constraints on the shape 
and normalization of the C iv column density distribu- 
tions at z = 5.3 — 6.0. For a power-law distribution, 
f{Nciv) oc N^"^^, with a slope equal to that measured 
at z ~ 3 (a « 1.8), the normalization must decline by 
at least a factor of 4.4 (95%). The distribution may also 
fiatten or steepen, but only a narrow range of normaliza- 
tions are allowed for each a by our data and the fact that 
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at least two strong (logiVciv > 13.8) sys tems are seen 
in the two previously s tudied sightlines (|Siincoel 120061 : 
iRvan- Weber et aI1l2006l and confirmed herein). Our lim- 
its on fie IV are less strong due to the fact that the co- 
moving density can be domiated by rare, strong systems. 
Within the range 13.0 < logTVciv < 15.0, however, the 
constraints on f{Nc iv) imply that fie iv must decline by 
at least a factor of two from z ~ 2 — 4.5 to z = 5.3 — 6.0 
(see also Ryan- Weber et al., in prep). The fact that 
f{Nc iv) and fie iv are relatively constant at z < 4.5 fur- 
ther means that most of the buildup of C iv in the IGM 
must occur over an interval Az < 1, or At < 300 Myr. 

The fact that C iv in the IGM rapidly increases be- 
tween z ^ 5.3 and z ~ 4.5 suggests that ionization ef- 
fects may be playing an important role. That is, a sig- 
nificant fraction of intergalactic metals we see as C iv 
at z 4.5 may already be in place at z > 5.3, but in a 
lower ionization state. This may especially be true if the 
metals are confined to relat ively high overdensities (e.g., 
iQppenheimer fc Davell2006l ). but other factors may also 
be at work. An increase in the overall IGM metallicity, 
a wider dispersal of metals by galactic winds, and an 
increase in the C iv/C ratio due to an increase in the 
intensity or a hardening of the UV background would all 
help to explain the rapid increase in C iv at z < 5.3. 
To evaluate the relative importance of these effects, it 
is critical to measure the abundance of metals in lower 
ionization states. We will do this in a subsequent work. 

The contrast between our lack of detections along 
four sightlines and the presence of two strong C iv 
absorbers towards SDSS J1030-I-0524 suggests that 
there may be large-scale variations in IGM enrichment 
and/or in the ionization state of the metal-enriched 
gas. We are currently planning follow-up observations 
of SDSS J1030-I-0524 to look for weak C iv systems. 
The presence of additional C iv would more strongly 
indicate that this sightline is either more enriched or 



more highly ionized than those in our current sample. 
Alternatively, the strong C IV systems may be associ- 
ated with rare, UV-bright star-forming galaxies, in which 
case the enhancement of C iv would be a local effect. 
This is supported by the overabundan ce of z-dropouts 
in the fiel d of SDSS J1030-H0524 (Sti avelh et all 120051 : 
iKim et aLl 20081. although spectroscopic follow-up is re- 
quired to be certain that these galaxies are associate with 
the C IV systems. 

Finally, we note that the decline of C iv at z > 5.3 
poses an interesting challenge to the use of low- ionization 
metal lines as tracers of reionization. If the number den- 
sity of C IV systems had remained constant out to z ~ 6, 
then we would be more certain that the IGM is suffi- 
ciently enriched at early times for metal lines to serve as 
meaningful probes of the ionization state of the gas. The 
disappearance of C iv at high rcdshifts, however, means 
that if an "O i forest" is not found at z > 6, the result 
could be attributed to a deficit of metals rather than to 
a highly-ionized IGM. On the other hand, the detection 
of a significant number of neutral metal absorption sys- 
tems at z ~ 6 would have a double significance: that the 
IGM is at least partially metal-enriched by that redshift, 
and that pockets of the IGM were not permanently ion- 
ized until well after the early reionization epoch favored 
by CMB studies. This will be the focus of the following 
paper. 
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